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Comment  on:  An  intermolecular  potential  for 
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and 
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1  Definitions  and  Correction 

A  recent  paper  by  Cappelletti  et  aZ.[l]  (CVPHD)  makes  use  of  classical  tra¬ 
jectory  (CT)  effective  kinetic  cross  sections  on  an  optimized  potential  energy 
surface  to  predict  Z  for  nitrogen  gas  over  a  range  of  temperatures.  It  is  sug¬ 
gested  here  that  an  alternative  method  of  conversion  between  the  calculated 
cross  sections  and  the  rotational  relaxation  quantities  (tjv,  the  volume  viscosity, 
and  Z,  the  rotational  collision  number)  may  give  more  physically  meaningful 
results. 

First,  a  clarification  regarding  a  typographical  error  in  eq.  (3)  of  CVPHD. 
Z  is  defined  as:  rrot/rc  where  Trot  is  the  rotational  relaxation  time  and  rc  is  the 
collision  time. 
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_  V*  (1) 

T'01  -  q2p  ~  p3®(000i)  ’ 

where  q 2  is  defined,  as  kc^Jcl  k  is  Boltzmann’s  con¬ 
stant,  p  is  the  pressure,  T  vs  the  temperature,  ins  the 
average  thermal  velocity,  and  <5(0001)  is  the  effective 
cross  section  for  first-order  internal  energy  relaxation, 
while  <5(0010)  is  the  cross-section  for  translational 
energy  relaxation.  The  approximate  result  is  true  to 
first  order  where  we  have  used: 

m  j&g  _ _ ££-,  (2) 

n’  ~  3cj5S(0010)  56(0001) 

6(0010)  =^6(0001).  (3) 

The  superscript  (1)  indicates  a  first  order  result  and 
the  superscript  (2),  used  below,  a  second-order  result. 
We  use:  c.  is  the  isobaric  ideal  gas  heat  capacity,  cu  is  the 
isochoric  heat  capacity,  and  c;m  is  the  internal  isochonc 
heat  capacity 


_  t*7)  xj  n 
Tc  —  4  p  4pi5S(2000)  ’ 

where  n  is  the  shear  viscosity.  In  the  case  of  shear  vis¬ 
cosity,  higher-order  corrections  are  quite  small.  Com¬ 
bining  expressions  leads  to: 

_  -4q»  45(2000)  m 

Z“^?^7rS(0001)' 

This  expression,  not  equation  (3)  of 
duces  the  first  order  2  curve  of  figure  (6)  in  CVPHD. 
The  curve  in  that  figure  which  is  labelled  second-order  Z 
cannot  be  consistently  reproduced. 

The  expression  for  second-order  volume  viscosity, 
which  will  be  used  below  is: 

m  (6) 

^  “  3cjD  S0’ 

where  S°  (as  defined  by  Heck  et  al.  [2])  is  the  deter¬ 
minant  of  a  matrix  of  sixteeri  effective  cross-sections, 
but  is  most  strongly  influenced  by  <5(0010),  and  Sfr  is 
the  determinant  of  the  minor. 
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2  Comment 

The  point  of  this  comment  is  to  raise  the  question  of  the  most  appropriate  value 
of  dm  to  use  in  the  above  expressions.  Table  1  shows  the  values  for  cp,  the  iso- 
baric  ideal  gas  heat  capacity,  and  the  implied  values  for  c„,  the  isochonc  heat 
capacity,  and  dm,  the  internal  isochoric  heat  capacity,  as  well  as  the  values  for 
shear  viscosity,  recommended  by  Millat  and  Wakeham  [3]  for  several  tempera¬ 
tures.  It  has  been  assumed  by  CVPHD  that  the  values  for  dm  due  to  Millat 
and  Wakeham  should  be  used  in  the  above  expressions.  If  this  is  done,  then 
the  identity  eq.  3  does  not  hold  true,  and  one  may  obtain  two  different  results 
for  tjv  and  Z  by  using  the  CT  result  for  6(0010)  or  by  using  the  CT  result  for 
6(0001)  and  eq.  3  to  calculate  6(0010).  Row  A  of  the  table  shows  the  results 
when  6(0001)  is  used  with  eq.  3  to  obtain  a  value  for  6(0010)  which  is  then 
put  into  eqs.  2  through  6  to  obtain  r)v  and  Z  in  first  and  second  order.  Row  B 
shows  the  result  when  6(0010)  is  used  directly  in  eqs.  2  through  6. 

However,  for  temperatures  above  300  K,  the  measured  value  of  dnt  begins 
to  reflect  an  appreciable  contribution  from  the  vibrational  mode  of  the  nitrogen 
molecule.  This  contribution  is,  of  course,  not  accounted  for  in  the  calculated 
cross  sections,  which  assume  a  rigid-rotor  molecule.  Therefore,  it  is  inconsistent 
to  combine  a  dm  =  d-ot+Cvib  value  in  equations  with  cross  sections  that  are  for 
rotation  only.  It  is  preferable  to  use  dm  =  1  in  the  above  expressions,  which  is 
the  value  that  theoretically  would  be  observed  (above  300  K)  if  the  vibrational 
mode  were  rendered  completely  inactive,  to  be  consistent  with  the  rigid-rotor 
assumption  of  the  CT  calculations. 

Of  course,  if  one  does  this,  one  should  be  careful  not  to  compare  the  value 
obtained  for  volume  viscosity  with  any  measured  value  that  includes  a  significant 
contribution  from  the  vibrational  mode.  Comparisons  should  still  be  valid  with 
any  measured  quantities,  or  model  predictions,  that  include  only  contributions 
from  rotational  relaxation.  This  would  include  predicted  values  for  Z  based 
on  state-to-state  rotational  energy  transfer  models  [4];  or  measured  values  for 
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volume  viscosity  if  the  rotational 
relaxation  time  is  so  much  shorter  than  the  vibrational 
relaxation  time  that  it  can  be- measured  as  a  separate 
quantity.  Such  comparisons  would  still  require  the 
assumption  that  the  rotational  mode  is  completely 
decoupled  from  the  vibrational  mode.  Note  that  Carne- 
valc  ei  al.  [5]  specifically  assume  that  vibration  is  frozen  in 
their  analysis  of  nitrogen  acoustic  absorption  data,  even 
up  to  their  maximum  temperature  of  1300K. 

The  results  obtained  when  the  value  c^,  =  crel  =  J i 
and  cv  =  2.5k  are  used  in  the  above  expressions  are 
shown  in  row  C  of  the  table.  In  this  case,  no  additional 
column  is  needed  to  show  the  results  based  on  6(0001) 
or  S(001 0),  since  for  cjn,  =k  the  identity  equation  (3) 
holds  true  and  both  cross-sections  yield  the  same  values 
of  77„  and  Z.  c„  =  2.51c  is  the  value  that  would  hypo¬ 
thetically  be  observed  if  vibration  were  rendered  com¬ 
pletely  inactive,  which  is  the  assumption  most  consistent 
with  the  rest  of  the  present  argument.  Note  that  the 
results  in  row  A  are  rather  fortuitously  close  to  those 
in  row  C.  This  is  because  the  use  of  the  measured  value 
for  Cjni  increases  the  term  c?n{/c*  in  equation  (2),  but  the 
measured  value  for  then  produces  an  offsetting 
decrease  due  to  a  larger  value  of  S(0010)  based  on  the 
relation  in  equation  (3). 


I 

i  Figure  1.  Second-order  results  for  2  using 
j  the  three  methods  of  calculation:  Case 
A,  (5(0001  )CT;  Case  B, 
j  ©(OO10)CT;  Case  C,  -  -  -  ,  a  * 


(  The  three  different  results  for  Z  (second-order)  for  the 
’  temperatures  discussed  here  are  shown  in  figure  1.  The 
.  result  for  row  B  can  be  seen  to  turn  up  sharply  at  higher 
|  temperatures.  This  behaviour  would  appear  to  be  incon- 
i  aistem  with  the  best  informaiion  currently  available  on 
1 Z  for  nitrogen  from  a  variety  of  sources  [6].  In  addition, 
the  use  of  the  measured  value  of  creates  the  difficulty 
of  arbitrarily  choosing  one  of  the  two  calculated  cross- 
i  sections  which  produce  conflicting  values  of  Z.  These 
.observations  support  the  consistency  argument  put  for¬ 
ward  above  for  using  c,nL  =  crol  =  k  in  this  context.  In 
! addition,  the  need  to  move  toward  a  method  of  calcu¬ 
lating  transport  cross-sections  that  is  not  limited  to  the 
rigid  rotor  model  is  underlined. 

*  v 
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H 


temperature 

300  K 

500  K 

1000  K 

1500  K 

2000  K 

measured 

Cp/fc 

3.507 

3.560 

3.934 

4.192 

4.328 

heat 

Cv/k 

2.507 

2.560 

2.934 

3.192 

3.328 

capacity 

Cint/k 

1.007 

1.060 

1.434 

1.692 

1.828 

shear  viscosity 

V 

17.9 

26.0 

41.2 

53.9 

65.7 

Case  A 

rjv 1} 

9.7 

17.5 

35.9 

51.4 

65.7 

11.3 

20.9 

41.9 

59.0 

74.8 

6(0010)  = 

z w 

4.3 

5.3 

6.7 

7.3 

7.7 

2Cr0t/k& (0001)  CT 

^(2) 

5.0 

6.3 

7.8 

8.4 

8.8 

Case  B 

9.8 

18.5 

51.4 

87.0 

120.2 

Vv2) 

11.4 

22.4 

64.6 

111.1 

155.1 

6(0010)  = 

ZW 

.4,4 

5.6 

9.5 

12.4 

14.1 

6(0010)ct 

Z&) 

5.0 

6.8 

12.0 

15.8 

18.2 

Case  C 

9.8 

17.4 

34.3 

49.5 

63.8 

rP 

11.3 

20.9 

43.2 

63.3 

82.3 

Cint  —  1 

z^y 

4.4 

5.3 

6.6 

7.3 

7.7 

P 

111 

to 

01 

z{2) 

5.0 

6.4 

8.3 

9.3 

9.9 

Table  1:  Calculated  values  for  t]v  and  Z  using  CT  cross  sections.  Values  for  77 
and  r]v  are  in  units  of  10-5<7/(crosec). 
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